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Scientific community has well recognized that a Lorentzian medium exhibits anomalous dispersion 
behavior in its resonance absorption region. To satisfy the Krammers-Kronig relation, such an anomalous 
region has to be accompanied with significant loss, and thus, experimental observations of negative group 
velocity in this region generally require a gain-assisted approach. In this letter, we demonstrate that the 
negative group velocity can also be observed in the absence of absorption resonance. We show that the 
k-surface of a passive uniaxial Lorentzian medium undergoes a distortion near the plasma frequency. This 
process yields an anomalous dispersion bandwidth that is far away from the absorption resonance region, 
and enables the observation of negative group velocity at the plasma frequency band. Introducing 
anomalous dispersion in a well- controlled manner would greatly benefit the research of ultrafast photonics 
and find potential applications in optical delay lines, optical data storage and devices for quantum 
information processing. 

n electromagnetics, the classical Drude-Lorentz model offers an intuitive and accurate description of both 
normal and anomalous material dispersion in an extended spectrum from direct current to optical frequen- 
cies 1,2 . As a direct consequence of the Kramers -Kronig relation, anomalous dispersion occurs as the 
wavelength approaches to the electric scale of the absorption resonance of the medium 3 , exhibiting a reversed 
spatial-temporal frequency dependence in the opposite of the normal dispersion case. Without violating the 
causality 4 7 , anomalous dispersion in the absorption resonance region implies an electromagnetic (EM) wave 
packet propagating with a superluminal or negative group velocity, i.e., an EM pulse appears to exit the medium 
before it enters 817 . In this paper, instead, we demonstrate that the negative group velocity can also be observed in 
the absence of absorption resonance. We show that when an EM wave is incident obliquely upon the surface of a 
uniaxial medium with a Lorentzian dispersion along its optical axis, the k-surface of the refractive wave undergoes 
a distortion as the frequency sweeps across the plasma frequency. Such a process yields an anomalous dispersion 
bandwidth that is far away from the absorption resonance region, and enables a new scheme for dispersion 
engineering in a well-controlled manner. This paper is structured as follows. We will first show the phase 
matching analysis of the proposed approach; next, we will report on our full wave simulation results. We will 
then report on our experimental observation of the negative group velocity at the plasma frequency band. Finally, 
we will draw our conclusion at the end of this paper. 

Results 

Phase matching analysis. The group velocity describing the propagation of wave packets in an isotropic medium 
can be expressed as v s = c/n g , where n g = n + co-dn / dco, n is the usual refractive index, c is the speed of light in 
vacuum and CO is the angular frequency. In the anomalous dispersion region with strong absorption resonance, 
the term dn/dw in n g is negative, and as such, the group velocity can be a negative value if this term is sufficiently 
small. This is exactly the phenomenon that has been reported previously 4-8-17 . Normally, this phenomenon does 
not occur at the plasma frequency band, which is a normal-dispersive region in the Drude-Lorentz model. 
However, the situation is different when the medium is anisotropic 1819 . Recently, an optical topological 
transition in a strongly anisotropic metamaterials has been described 20,21 . In this paper, we further point out 
that during such topological transition, the k-surface undergoes a distortion. This process yields an anomalous 
dispersion bandwidth and enables the observation of negative group velocity at the plasma frequency band. 

Without losing generality, we begin with a simplest uniaxial medium with a scalar permittivity e r = 1 and a 
permeability tensor Ji r = diag(l > l,lA n ), where ji rz = 1 — CQp/[co(co + iy e )] follows a frequency dependence of Drude 
model with a plasma frequency w p and a damping frequency y e . Notice that all of above constitutive components 
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are normal dispersive and, unlike the absorption resonance case, the 
loss of j.i rz characterized by y e can be very small, as indicated later. 

As shown in Fig. la, with an oblique EM incidence from the air 
(electric field polarized along the y direction) and a wave vector k = 
xk x + zk z , the dispersion relation (or k-surface) in the medium is 

k 2 Jti rz + kl = i, r kl = RAm/c) 2 , (1) 

where k 0 is the wave number in air, k x and k z are wave numbers in 
x and z directions, respectively. Based on the phase matching 
condition, for a given incident angle 0 , we have k x = k 0 sinO and 

k z = + \] £ r (a>/c) 2 — k\l pi rz . Due to the passivity nature of the dis- 
cussed medium, the sign of k z can be determined by choosing the k z 
with a positive imaginary part. 

For a specified damping frequency, the shape of the k-surface 
described by equation (1) would strongly depend on the sign of 
fi rz , as shown in Fig. lb. In a hypothetical lossless case shown in 
the upper panel of Fig. lb (the damping frequency is set to 0), the 
k-surface could be either a closed ellipse (cu > oj p , fi rz > 0) or an 
open, two-branch hyperbola (w > w p , fi rz < 0), and at the exact 
plasma frequency with ^ rz = 0, the ellipse and the hyperbola shrink 
separately to a line segment (-Fi-p2> when w —»<»„+) and two rays 
(FiRi and F2R2, when CO —> to p —) on the z axis. Therefore, under the 
lossless hypothesis, there is a discontinuity of the k-surface at the 
plasma frequency: when the frequency sweeps across the plasma 
frequency, the k-surface changes abruptly from a line segment to 
two rays, which breaks the continuity of dispersion and implies the 
non-causality of the lossless hypothesis 22 . On the other hand, when 
the damping frequency is not zero (shown in the lower panel of 
Fig. lb), the k-surface distorts from a closed ellipse to an open, 
two-branch curve that looks like a hybrid of an ellipse and a hyper- 
bola. The continuous evolving of the k-surface from a distorted 



ellipse to a distorted hyperbola implies the causality and physical 
reliability of such a Lorentzian medium. 

Next we focus our discussion on the continuous distortion of the 
k-surface at the vicinity of the plasma frequency. In Fig. lc, we show 
the phase matching diagram for two adjacent frequencies (co and to 
+ 8) in the transition zone (i.e. around the plasma frequency) of the 
k-surface distortion. We define a critical angle 0 C (see Supplementary 
Note 1) for the discussion of the following two cases: 1) 61 < 9 a an 
inclined EM wave propagates at the forward direction with a group 
velocity of v g i; 2) 0 2 > 0 C , an inclined EM wave propagates at the 
backward direction with a group velocity of v g2 - We also note that, 
the phase velocities denoted by v pi andv^, 2 remain positive, meaning 
that the phase velocities of the carriers always propagate forwardly in 
both cases. Therefore, the negative group velocity v g 2 is not accom- 
panied with a negative refraction. From equation ( 1 ) , both v g \ and v S 2 
can be calculated by 



v g = Vj-CO(fc) = xv gx + zv„ z — x 



1 



1 



Re [dk x /dio] Re [dk z /doj] ' 

In the rest of the paper, we focus our discussions on the group 
velocity at the propagating direction, i.e. the z direction. As shown in 
Fig. Id, for a smaller incident angle 0 t (green curves), k z is normal- 
dispersive with a small loss term denoted by the dash green curve; 
while, for a larger angle 0 2 , k z exhibits an anomalous dispersion 
behavior inside the transition zone (marked in grey, also see 
Supplementary Note 2). We note that, regardless the normal-dis- 
persive property of j.i rz (the inset of Fig. Id), the dispersion of k z is 
anomalous, and the behavior of k z in this transition zone is similar to 
that in the conventional absorption resonance region. 

Full wave simulation. We performed a full wave EM simulation 
based on the time domain propagation of wave packets using a 
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Figure 1 | Theoretical analysis, (a) Diagram of EM wave packets impinge on the surface of a uniaxial plasma medium, (b) Illustration of the k surfaces 
around plasma frequency in cases with (bottom) and without loss (top), (c) Phase matching at the boundary between air and the discussed medium. 



(d) Dispersion of k z for different incident angles, where k p - 
0.02co„. 



pic, Inset shows the dispersion of /i rz in the same frequency zone, where y e is assumed to be 
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Figure 2 | Full wave simulation of group velocities in a uniaxial plasma medium with oj p = 10 GHz and y e = 0.02ra p for different carrier frequencies. 

(a) The frequency of carrier is 10 GHz and (b) The frequency of carrier is 12 GHz, where the incident angle is 10°. The blue curve denotes the detected 
incident EM packets at the surface of the discussed medium. The amplitude of each waveform is normalized by the incident wave. 



commercial software package, CST microwave studio. As shown in 
Fig. la, we recorded the signals of the wave packets in a 20-cm-thick 
magnetically uniaxial medium. The optical axis of this medium is at 
the z direction and the dispersion of \i rz follows the Drude model, 
with ojp = 10 GHz and y e = 0.02(O p . The incident angle was chosen 
to be 10 degrees, larger than the critical angle 0 C which is determined 
to be 1.8 degrees under this configuration. In order to measure the 
signal amplitudes inside the medium, three probes were placed at 
point 'a', 'b' and 'c', with a 70-mm increment, as shown in Fig. la. 

In this case, the transition zone of this configuration is estimated to 
be 9.93-10.22 GHz (see Supplementary Note 2). Two simulations 
were performed for frequencies both inside (Fig. 2a) and outside 
(Fig. 2b) the transition zone. In the first simulation, a 10-GHz carrier 
modulated with a 40-MHz sinusoidal envelope signal was simulated 
as the incident wave. The detected signals at the three points, as well 
as the incident wave at the surface of the medium, are shown in 
Fig. 2a. Several distinct features are worth noting: 1) due to the 
imaginary part of k z in Fig. Id, the signal intensity decays with an 
attenuation rate of — 1.28 dB/cm from point 'a' to point 'c'; 2) for the 
first wave packet, the 'front velocity' of the envelope signal is pos- 
itive 23 ; 3) except the front knot of the first packet, all the successive 
peaks and knots are first seen at the furthest point 'c' (see the mag- 
nified panels labeled with T' and '2'), implying that they propagate 



backwardly from point V to point 'a'; 4) Compared with the incident 
wave measured at the surface (blue line), the wave front of the second 
wave packet arrives at point 'c' even before the first wave packet 
entirely enters into the medium, implying the existence of the nega- 
tive group velocity. We note that, despite the fact of the negative 
group velocity, the front velocity of the signal is always positive, in 
other words, the information velocity carried by the wave front is still 
positive. As a comparison, we performed the second simulation with 
a carrier frequency of 12 GHz in Fig. 2b as a reference. In this case, 
the frequency is not in the transition zone (i.e. normal-dispersive 
behavior is expected). It is seen that the front and group velocities 
are all positive in Fig. 2b. 

Experimental verification. We demonstrated the experimental 
verification of our analysis based on the concept of metamaterial 24 . 
Consisting of densely arranged artificial, sub -wavelength electric 
and/or magnetic dipole moments, metamaterials are capable of 
providing predefined effective constitutive parameters based on the 
Maxwell-Garnett effective medium theory 25,26 . Fig. 3a shows the slab- 
like metamaterial sample. Consisting of multilayer printed circuit 
boards carrying periodically arranged sub-wavelength split-ring 
resonators (SRRs), the dispersion of the effective z-component 
permeability (fi rz , shown in the lower inset) closely follows a 
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Figure 3 | Experimental structure and setup, (a) Metamaterial sample and the real and imaginary parts of ji rz retrieved from simulation data, 
(b) Experimental setup. 
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typical Drude-Lorentz dispersion model with an absorption 
resonance around 9.19 GHz and a typical plasma dispersion 
region between 9.35 GHz and 9.65 GHz (marked in gray, with 
ojp~ 9.45 GHz and a small imaginary part varying from 0.02 to 
0.2). Fig. 3b shows the experimental setup for generating and 
detecting the EM wave packets. Detailed information of the 
metamaterial sample and the experimental setup are provided in 
Methods. 

With a coherence measurement of the transmitting wave packets 
inside the sample, the demodulated envelope signals with the carrier 
frequency of 9.5 GHz (in the transition zone) and 10 GHz (out of the 
transition zone) were directly obtained, as shown in Fig. 4a and 4b. In 



Fig. 4a, the amplitude of envelope decays as the wave propagates from 
point 'a' to point 'c' with an attenuation rate of —2.18 dB/cm, cor- 
responding to a propagation length of 1.99 cm. On the other hand, 
the knots (as well as the peaks) of wave envelop are first seen at the 
furthest point 'c', and they propagate backwardly to point 'b' (1.51 ns 
later) and then point 'a' (2.57 ns later). As a comparison, in the 
10-GHz carrier case (outside the transition zone), the knots of wave 
envelop propagate from nearest point 'a' to point 'b' (0.83 ns later) 
and then point 'c' (1.58 ns later) with an attenuation rate of 
—0.16 dB/cm, corresponding to a propagation length of 27 cm. 
Fig. 4c shows the measured group delay Af (upper panel), the attenu- 
ation rate (lower panel) versus carrier frequency in the transition zone 
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Figure 4 | Experimental results and analysis, (a) The demodulated envelope signals in time domain measured at points 'a', 'b' and 'c', where the 
frequency of the carrier is 9.5 GHz. (b) The same measurement when the frequency of carrier is 10 GHz. (c) The measured and theoretical group delay At 
(upper panel) and attenuation (lower panel) versus carrier frequencies in the transition zone with an incident angle of 15 degrees, (d) The measured and 
theoretical group velocity v— versus carrier frequencies, (e) The measured and theoretical group velocity v gz versus different incident angles at 9.5 GHz. 
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with an incident angle of 15 degrees. For comparison, we also show in 
Fig. 4c the group delay and attenuation rate calculated by equation (2) 
with the simulation-based e r and ji rz shown in Fig. 3a. We see that 
both experimental results of group delay and attenuation rate are 
consistent with the trend of the theoretical calculations, except for a 
slight difference in central frequency and bandwidth. Comparing 
with the experimental results, the calculated bandwidth is 40 MHz 
narrower with a frequency shift of 30 MHz. This slight discrepancy 
can be attributed to the difference between the actual and simulation- 
based e r and fi rz . It is well known that the loss of a fabricated meta- 
material is higher than the simulation case due to the misalignment of 
unit cells and the incontrollable machining errors (In Supplementary 
Note 3, it is shown that either a higher material loss or a larger 
incident angle leads to the increase of the bandwidth of the transition 
zone). 

In Fig. 4d, the z component of the group velocity is calculated by 
Vgz = Ad / At , where Ad is the distance between point 'b' and point 'c', 
i.e., 48 mm. Notice that the group velocities marked in blue dots, 
dashed red line and solid green line are calculated from the corres- 
ponding points and curves in Fig. 4c. As a direct consequence, v gz is 
negative in the transition band. The calculated curves also show that 
v gz approaches to negative infinity at the edges of the transition zone, 
reflecting the fact that the group delay At continues at zero when v gz 
changes its sign. Fig. 4e shows the dependence of the group velocity 
Vgg versus incident angle at the carrier frequency of 9.5 GHz. We see 
that v gz is positive at small incident angles (0-5 degrees) and negative 
at larger incident angles. A negative group velocity as fast as —0.3c is 
measured at 9.5 GHz with an incident angle of 7.5 degrees. 

Discussion 

From Fig. 2, we see that the material loss, the anisotropy and the 
oblique incidence play indispensable roles in the generation of the k- 
surface distortion, which finally results in the anomalous dispersion. 
This particular k-surface distortion cannot appear in an isotropic 
medium, nor in the case with a normal incidence, which makes 
our approach differentiate itself from previously reported ones. 
However, causality still holds: in Fig. 3, both the front velocity and 
the group velocity (reflecting the phase velocity of the wave packets) 
are simultaneously illustrated, and we see that although the wave 
packets travel at a superluminal velocity, the information carried 
by them is not superluminally transmitted. In Fig. 4, both the fre- 
quency and the incidence angle dependence of the negative group 
velocity are demonstrated, which provides a new degree of freedom 
to engineer the dispersion property. It is also worth noting that, with 
an oblique incidence, the dispersion of group velocity would spread 
the wave packets not only in temporal axis, but also in spatial axes. At 
an incidence angle where the dispersion of group velocity is strong 
enough, e.g., at the critical angle 0 C , multiple carriers modulated with 
different information can be spatially separated, and as such, it 
provides a potential approach for a temporal-spatial information 
processing. 

To summarize, we have demonstrated the introducing of anoma- 
lous dispersion in a passive medium. We experimentally measured 
an unusual travelling of EM wave packets with a negative group 
velocity at the plasma frequency band, which is far away from the 
absorption resonance region. This abnormal wave packets propaga- 
tion arises from the physical requirement of the dispersion continu- 
ity at the plasma frequency. Different from previously reported 
superluminal experiments, our approach does not require strong 
absorption resonance, and can be directly measured inside the 
medium. The proposed approach for engineering group indices in 
normal dispersive media is considerably easier than the existing 
methods. Introducing anomalous dispersion in a well-controlled 
manner would also greatly benefit the research of ultrafast photonics 
and is therefore suitable for diverse applications such as optical delay 
lines, optical data storage and devices for quantum information 



processing. Recent progress of active nanoplasmonic metamtaerials 
shows that active methods can be used to effectively compensate the 
metal losses, providing the possibility for fabricating low-loss optical 
uniaxial metamateirals in the future 27,28 . 

Methods 

Fabrication of the experimental sample and the retrieval of [i rz . The slab-like solid- 
state metamaterial sample consists of multilayer printed circuit boards (PCBs) 
carrying sub-wavelength split-ring resonators (SRRs). In the fabrication process, the 
metallic (copper) SRRs were printed on 0.5-mm-thick FR4 substrates, with a relative 
permittivity around 4.6 at microwave frequencies. Each board contains 82 unit cells 
along the x direction with a periodicity of 6 mm and 15 unit cells in the y direction 
with a periodicity of 10 mm. Eleven layers of such PCBs were aligned along the z 
direction spaced by 16 mm to form a slab-shaped sample with a width of 500 mm 
along the x direction, 150 mm in the y direction and a thickness of 160 mm in the z 
direction. 

For an oblique incident wave with the electric field polarized in the y direction 
impinging upon such a medium, magnetic resonances can be induced at the metallic 
rings, yielding a strong dispersion of permeability with a negative band of fi rz . 
Meanwhile, there is no magnetic response in x and y directions, and thus /t rA - and ^ ry 
can be treated as unity. By utilizing the homogenization approach 29,30 , the dispersion 
of /i rz can be retrieved from the scattering parameters from a full-wave simulation, 
and the result is shown in the inset of Fig. 3a. We see that the dispersion of pi rz closely 
follows a classical Drude-Lorentz dispersion model: there is a strong absorption 
resonance around 9.19 GHz, and there is a typical plasma dispersion region between 
9.35 GHz and 9.65 GHz (marked in gray) with w p ~ 9.45 GHz and y e ~ 0.05co p . 
Different from the previous unity assumption, the retrieved s r is around 2. Utilizing 
the retrieved \i rz and e r , the transition zone of the k-surface distortion can be estimated 
as from 9.42 GHz to 9.51 GHz with a 15-degree incident angle. We note that, in the 
transition zone, the pi rz is normal-dispersive and its imaginary part is quite small 
(varying from 0.02 to 0.2). As a comparison, the imaginary part is up to 7 in the 
conventional absorption resonance region). 

Experimental setup. The experimental setup is shown in Fig. 3b. A 10 MHz 
sinusoidal envelope signal was generated by a signal source (Agilent's 33210A) and 
used to modulate a carrier provided by a microwave source (Agilent's E8257D). The 
modulated signal was amplified and sent to a horn antenna located 150 cm away from 
the sample. The EM wave leaving from the antenna was electrically polarized along 
the y direction, and impinged on the surface of the sample with an incident angle of 
15 degrees. We placed the probe (a monopole antenna) at points 'a', 'b' and 'c', 
separately, to measure the intensity of the wave packets inside the sample. The 
detected signals from the probe were then demodulated by a crystal detector (Yamei 
Microwave's BO-20A) to recover the envelope signal. Triggered by the primitive 
envelope signal generated by the Agilent 33210A, the recovered envelope signal was 
displayed and sampled using a digital oscilloscope (Tektronix's MSO 4104B). To 
remove high frequency noise and waveform jitter, the built-in 20-MHz low pass filter 
and the average function of the MSO 4104B were enabled to obtain steady and clear 
envelope waveforms (see Supplementary Fig. S3 and S4). 
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